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1.   Introduction 

Experimental  knowledge  of  the  flow  field  generated  by 
rotating  turboimpellers  is  essential  for  the  research  and 
development  of  turbomachinery .   This  information  is  used  to 
refine  design  methods,  develop  new  flow  models  which  include 
secondary  flow  and  tip  clearance  effects,  and  especially  to 
verify  computer  programs  designed  to  calculate  flow  throuqh 
rotating  blade  rows. 

Laser  velocimeters  have  been  used  successfully  in  recent 
years  to  measure  the  flow  inside  and  downstream  of  rotors 
(see  Ref .  1) .   Certain  disadvantages  have  become  apparent, 
however.   The  laser  techniques  are  reliable  only  in  the  hands 
of  experienced  investigators,  the  pressure  field  remains 
unknown,  and  usually  the  measurement  of  more  than  two  compo- 
nents of  the  velocity  field  is  complicated  and  expensive. 
Furthermore,  it  is  difficult  to  perform  measurements  close  to 
walls.   Development  of  alternative  techniques  to  overcome 
these  deficiencies,  as  well  as  to  achieve  redundancy  in  mea- 
suring the  flow  field,  are  reasonable  and  worthwhile  tasks. 

This  report  describes  a  particular  method  and  the  com- 
putational support  necessary  to  measure  the  flow  field  behind 
an  impeller  in  the  stationary,  bladeless  gap. 


2.   Description  of  Method 

The  following  method  requires  two  semiconductor  pressure 
probes  along  with  a  technique  for  synchronized  sampling  for 
determining  the  fluid  velocity  vector  downstream  of  a  rotor. 

The  two  probes  (see  Fig.  1)  are  positioned  inside  the 
machine  casing  so  they  will,  in  turn,  intercept  periodically 
the  same  part  of  the  flow  leaving  a  particular  passing  rotor 
passage.   Each  probe  reading  is  sampled  when  the  designated 
blade  passage  reaches  a  desired  position  relative  to  the 
probe.   Synchronization  is  achieved  through  a  suitable  method 
(Ref.  2,3). 

Four  quantities  are  needed  to  determine  the  velocity 
vector:   yaw  angle,  pitch  angle,  static  pressure  and  total 
pressure.   Accordingly,  four  measurements  must  be  made  to 
evaluate  these  unknowns.   By  rotating  the  probes  about  their 
tips,  pressure  readings  in  four  different  directions  can  be 
taken,  and  the  data  used  to  calculate  the  velocity  vector. 
Computer  program  VELOCITY,  given  in  Appendix  II,  was  developed 
to  perform  the  somewhat  arduous  calculations. 

The  geometries  of  the  two  probes  are  shown  in  Fig.  1. 
Before  being  used,  the  probes  must  be  calibrated  so  their 
responses  to  flows  coming  from  different  directions  are  known. 
A  highly  directional  probe  is  desired  to  increase  the  accuracy 
in  finding  the  yaw  and  pitch  angles,  and  consequently  the 
velocity  magnitude.   The  following  method  is  recommended  for 
calibrating  each  probe  - 
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Figure  1.   A  type  and  B  type  probes 


1.  Establish  a  steady,  controlled  flow  of  fluid,  and 
determine  the  velocity  vector  at  a  certain  region  of 
the  flow. 

2.  Position  a  probe  in  the  flow  and  rotate  the  tip  so 
that  a  sequence  of  pressure  readings  are  taken  for 
a  constant  yaw  angle  and  a  varying  pitch  angle. 
Repeat  the  procedure  at  a  new  yaw  angle  using  the 
same  pitch  angles.   The  result  will  be  an  array  of 
pressure  readings  corresponding  to  a  set  grid  of 
yaw  and  pitch  angles  (Fig.  2) . 

3.  From  the  known  flow  velocity  and  pressure  readings, 
a  coefficient  of  pressure  can  be  calculated  for  each 
angle  set: 


p-p  C   =  Coefficient  of  pressure 

=  where :  " 

P        Pm~P  ^ • 

c        rip  irs  p  _   pressure    reading 


p   =  static  pressure  of  flow 
p_  =  total  pressure  of  flow 

The  table  of   C  's   as  well  as  the  yaw  and  pitch  angles  which 
correspond  to  them  are  now  in  the  form  required  for  input  to 
program  VELOCITY. 

The  probe  calibrations  should  be  insensitive  to  Mach 
number  and  pressure,  and  are  not  valid  for  supersonic  flows. 
Should  any  significant  variations  in   C   be  observed  for 
different  flow  conditions,  further  calibrations  will  be  re- 
quired and  an  additional  iteration  scheme  added  to  the  com- 
puter program. 
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Figure  2.   Grid  of  Yaw  and  Pitch  Angles 


Experience   with    the   two-probe    technique   has    shown    that 
excellent   results    are   achieved  when    a  probe   type   A   is    rotated 
to    the    three   positions   +25°,    0°,    -25°    yaw. at   0°    pitch) ,and 
probe    type   B    is    used   at   0°    yaw   and   25      pitch,    Fig.     3) . 

The    two-probe    technique    is    strictly   applicable   only    to 
periodic    flows.      However,    data   obtained  on   successive    rotations 
of    the    rotor   can   be   averaged   to   eliminate   non-periodic    fluctua- 
tions.     This   was   effective    for    tests    reported   in    Ref.    2.,    where 
a   single   probe  was    used   to   establish    the   peripheral   blade-to- 
blade    distribution   of    flow  yaw   angle. 

It   is   noted   that   the   method   reported  here    is    a    further 
development  of    that   reported  earlier   in    Ref.    6,    and   overcomes 
some   of    the    earlier   limitations. 


PITCH 


DJ 


■Q 65 — ~ G — **  YAW 


O  PROBE  TYPE  A 
D  PROSE  TYPE  B 


•'i 


Figure  3.  Orientation  angles  of  the  probes  relative 
to  the  laboratory 


3.   Theory 

The  velocity  vector  for  a  three-dimensional  flow  can  be 
described  with  three  scalar  quantities.   The  nature  of  the 
problem  suggests  using  two  angles  (a  yaw  angle  and  a  pitch 
angle) ,  and  the  magnitude  of  the  velocity  (Fig.  4) . 

Since  pressures  and  not  the  velocity  are  measured,  the 
static  and  total  pressures  must  first  be  determined,  and 
Eq.  (1)  used  to  evaluate  the  velocity. 


El-  (i  +  IzImV^"1  (1) 

Ps         2 


Altogether,  four  unknowns  need  to  be  evaluated:   the  yaw  and 
pitch  angles,  and  the  total  and  static  pressures. 

Four  equations  are  needed  to  determine  the  four  unknowns, 
They  are  derived  from  the  four  pressure  readings,  each  pres- 
sure reading  having  been  taken  in  a  different  direction  as 
described  above.   The  following  equations  for  the  coefficient 
of  pressure  can  be  written: 


prps 

C  .  =  -= — -      i  =  1..4  (2) 

pi    PT"PS 


The  C  .'s  are  a  function  of  the  orientation  of  the  probe 
pi  F 

relative  to  the  flow;  i.e.  ,  for  a  given  flow  the  measured  C  ' 

P 

will   vary   measureably   as    the   probe    is    turned   into   and   away 
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Figure  4.   Velocity  Vector  V 


from  the  flow.  Each  "probe"*  will  have  its  own  C  charac- 
teristics  determined  experimentally.  The  result  will  be  a 
table  of  C   vs.  yaw  and  pitch  angles  for  each  probe. 

C  .  =  function  (an.  ,  <f>  . )       i  =  1..4      (3) 
pi  Ri    R  1 

For  realistic  problems,  only  one  point  (a  ,  <$>)    exists 

where  the  C  .'s  in  Eq.  (2)  will  equal  the  C  .,   of  Eq.  (3)  for 
pi       ^  M         pi's     ^ 

the  four  probes'  pressure  readings. 

The  probes'  characteristics  (C  's)  are  in  tabular  form 
because  they  cannot  be  represented  analytically  due  to  the  sten 
effect  and  production  inaccuracies.   Therefore,  a  numerical 
solution  to  the  problem  is  required.   The  procedure  chosen  for 
solving  the  problem  is  a  systematic  trial-and-error  search 
process,  essentially  a  convergence  scheme  on  two  variables: 
yaw  angle  and  pitch  angle. 

The  flow  direction  is  assumed  to  fall  within  some  set  of 
bounds,  defining  the  search  area  for  yaw  and  pitch  (Fig.  5). 
By  setting  up  a  grid  of  points  in  this  region  and  checking  how 
well  each  point  satisfies  the  criteria  of  equality  of  coeffi- 
cients of  pressure  (C  ,'s)  calculated  with  Eqs.  (2)  and  (3), 
the  point  with  the  smallest  error  can  be  found  and  used  as  a 
first  approximation  to  the  solution.   Repeating  this  procedure, 
only  with  a  smaller  grid  and  search  region,  will  result  in  a 
better  approximation.   This  sequence,  represented  in  Figs.  6 


*  Here,  the  term  "probe"  refers  to  a  particular  probe  type  in 
a  particular  position. 
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Figure  5.       Search  area 
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and  7  is  repeated  until  either  the  desired  accuracy  is  reached 
or  fatigue  sets  in.   Program  VELOCITY,  described  in  the  follow- 
ing section,  was  written  to  perform  these  calculations. 
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Figure  6.   Illustration  of  the  Search  Procedure 
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Figure    7.      Flow   Chart   of    the   Search   Procedure 
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4.   Program  VELOCITY 

Program  VELOCITY  was  written  to  perform  the  calculations 
outlined  in  the  previous  section.   A  description  of  the  program 
and  its  subroutines  is  given  below.   Fig.  8  summarizes  the 
major  sections  and  organization  of  the  program. 

For  each  run,  program  VELOCITY  reads  the  calibration 
tables  for  the  two  probes  from  files  outside  the  program. 
(Input  formatting  is  discussed  in  Appendix  V.)   Subroutine 
INPUT  performs  the  necessary  work,  and  can  be  modified  to 
accommodate  different  input  schemes  if  desired. 

The  fluid  temperature  and  molecular  weight  are  entered 
next.   These  properties  are  assumed  to  remain  constant  through- 
out the  run. 

The  settings  for  each  pressure  reading  are  read  next.   A 
setting  contains  the  following  data:   probe  type  (A  or  B) ,  yaw 
angle  setting,  and  pitch  angle  setting.   Again,  these  settings 
will  not  change  for  the  duration  of  the  run. 

Finally,  the  four  pressure  readings  are  entered. 

The  first  scan  is  initiated  and  covers  the  entire  region 
of  expected  flow  directions,  -40   to  +  40   in  both  yaw  and 
pitch  angles  in  the  present  case.   Points  are  chosen  every  5  , 
each  point  representing  a  unique  pair  of  yaw  and  pitch  angles. 
For  each  point,  a  static  pressure,  a  dynamic  pressure,  and  an 
error  are  calculated  by  the  scheme  described  below. 

A  point,  say  (a  ,  cj>)  is  tested;  i.e.,  a  test  is  performed 
to  prove  whether  assumed  flow,  oriented   a  degrees  yaw  and   $ 
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Figure  8.   Flow  Chart  of  Program  VELOCITY 
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degrees   pitch    relative    to    the    laboratory    reference    frame, 
corresponds    to    the    four  pressure    readings.      The    direction   of 
the    flow    relative    to   each   probe    setting   is    calculated.       For 
probe    setting      i    ,    oriented   at    (a.     ,    <J> . )    relative    to   the 
laboratory,    the   assumed    flow   approaches    at  a    relative    angle 
of: 

aRi    =   a-ou  (4) 


4>Ri  =    *-*±  (5) 


where  (a_..  ,  4>  . )  are  the  yaw  and  pitch  angles  respectively 

of  the  assumed  flow  relative  to  probe  setting   i  .   The  C 

P 

calibration  table  for  the  probe  used  in  setting   i   is  con- 
sulted and  a  C  (a„.  ,  <{>_.)  returned.   Subroutine  CPCAL  locates 

p   Rl      Kl 

or  calculates  the  desired  C  values  in  the  table.   The  scheme 

P 

used   in    CPCAL   is    a   search    technique    to    find    the   values   of   yaw 
and  pitch    surrounding   the    desired   point,    and    then    a   linear 
interpolation   over    these    four  points    as    shown    in    Fig.    9. 
Eq.     (2)    can   be    rewritten    in    the    form 

<Cpi)PT    +    d-Cpi^s    =    Pi  A    =    1"4  <6) 

the  only  unknowns  being  p    and   p   .   With  four  equations  and 
two  unknowns,  the  problem  will  be  inconsistent  unless  the  true 
a   and  $      were  chosen.   Accordingly,  the  following  schemes 
were  used  to  evaluate   p   ,  p    and  an  error. 
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Define: 


4 

C   =  E  C  (7) 

P    •  i      P- 


£  =  E  P,  (8) 

i=l  1 


C    =  minimum  of  (C    ,  C    ,  C    ,  C   ) 
pm  Pi    P2     P3     P4 


P__  =  p .  corresponding  to  the  C    chosen  above 
m     1         r  -o  p 

^m 


(C  )pT+  (4-C  )p   =  p  (9) 


and  also 


(C   )p_  +  (1-C   )p  =  p  (10) 

Prv,     T  P™     S         m 

^m  m 


These  two  equations  can  be  solved  for   pT   and  p   : 


p(l-C   )  -  p  (4-C  ) 
-    Pm     m   JE. 

pt  =  c  -  4c (11) 

P     Pm 


Ps  "      C  -  4cm  <12> 

P     P 
_£L      m 


4  p.-p 

Error  =  E   ABS(C i—^)  /4  (13) 

i=l       pi    Pt"PS 
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These  schemes  were  chosen  for  two  reasons: 

1)  They  used  all  the  available  data  to  derive  an  error 
which  would  effectively  represent  the  accuracy  of 
the  guess. 

2)  No  singularities  in  the  calculations  can  occur  except 
for  the  case  of  four  equal  C  's  (which  physically 
represents  trying  to  find  an  intersection  point 
among  four  parallel  lines) .   If  the  measurements  are 
taken  in  the  suggested  directions,  this  anomalous 
point  will  not  appear. 

For  each  point  guessed  in  the  initial  scan,  an  error  is 
calculated  and  the  point  with  the  smallest  error  is  saved.   A 
new,  finer  search  grid  is  composed  using  this  point  as  the  new 
origin.   The  boundaries  of  the  new  grid  are  the  points  from 
the  old  grid  which  were  closest  to  this  new  origin.   Referring 
to   Figure  10,  if   x   represents  the  true  solution,  the  new 
boundary  would  be  formed  by  the  points  marked  B-I,  and  the  new 
grid-width  would  be  one  third  as  large.   This  factor  was  chosen 
to  minimize  the  number  of  guess  evaluations.   (The  first  scan 
contains  a  large  number  of  guesses  in  order  to  correctly  isolate 
the  general  region  of  the  solution) . 

The  search  procedure  is  performed  on  each  new  grid,  and 
the  process  repeated  until  the  grid  width  is  less  that  0.5  . 
After  the  final  scan,  the  best  guess  is  used  to  calculate  the 
flow  velocity  and  Mach  number.   The  results  are  printed  out 
and  the  next  four  pressures  requested.   If  no  values  are 
entered  (end  of  data  set) ,  the  program  ends. 
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•*  TRUE  SOLUTION 

Figure  10.   Defining  new  grid  boundaries  from  the  nearest  neighbors 
of  the  point  with  the  smallest  error  "7 
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5 .   Discussion 

Extensive  tests  with  program  VELOCITY  have  led  to  the 
observations  and  suggestions  listed  below: 

1.  Excellent  results  are  achieved  when  the  probe  settings 
are  at  (yaw,  pitch)  angles  of  (-25,0),  (0,0),  (25,0) 
and  (0,25)  degrees.   This  corresponds  to  a  rotation 

of  probe  type  A  from  -25°  to  0   to  25  ,  and  one  read- 
ing from  probe  type  B  at  (0,25).   Poor  results  were 
achieved  for  the  symmetric  case  of  readings  at  (+25,0) 
and  0,+2  5)  degrees. 

2.  Highly  directional  probes  increase  the  accuracy  of 
the  procedure,  especially  if  the  C   variation  is 
significant  when  the  flow  is  nearly  head-on    To 
achieve  these  characteristics,  the  following  design 
suggestion  is  offered.   The  probe  can  be  formed  with 
a  spherical  tip,  the  pressure  tap  being  located  in 
the  center.   To  prevent  damage  to  the  sensitive  trans- 
ducer located  behind  the  pressure  tap  and  to  improve 
the  frequency  response,  the  void  between  the  pressure 
tap  face  and  the  transducer  should  be  filled  with  an 
an  appropriate  liquid  and  the  opening  of  the  pressure 
tap  sealed  with  a  thin,  low-inertia  membrane. 

3.  Higher  accuracy  naturally  results  if  more  calibration 
points  are  taken  for  the  probes'  C   tables.   The  linear 
interpolation  scheme  can  be  replaced  by  the  second 
order  scheme  offered  in  Appendix  5  (if  no  significant 
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anomalies   occur   in    the   calibrations) ,    the    second   order 
method    requiring    fewer    calibration   points    (say   every 
15    )     than    the    linear  method    (every    5      or   10    )  . 
4.       The    use    of   two   probes   of    relatively    simple    geometry 
in    periodic    flow   is    less    cumbersome    and    complex   than 
the   use   of    five-hole   probes    (Ref .    4) . 
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Notation  Summary 


C     -     Coefficient  of  pressure 


P 


Pi 


C   is  a  function  of   a   and  <$>    ,    C   =  C   (a  ,  <t>) 
P  P     P 


C  .    -     Coefficient  of  pressure  for  probe  setting   i 


C  .  =  C  (a_.  ,  $_.) 

pi     p   Ri  '  rRi' 


C  -  Sum  of   the    four   C    .'s 

P  pi 


C  -  Minimum  of    the    four   C    .    s 

P  pi 

^m  e 


Cpn      "  Cp(al    '    ♦l1 


P12  P1  2 


C  -  C    (a       ,    $-) 

p21  p      2  1 


C  -  C    (a       ,    4>    ) 

p22  p      2  2 


P  -  Pressure    (all   pressures    are    absolute) 

P.  -  Pressure    read    from  probe    setting   i 


p  -  Static   pressure 

pT  -  Total  pressure    (stagnation  pressure) 

p  -  Sum  of    the    four  pressures    (Pi's) 

P  -  Pressure    at   the    setting  where    C        occurred    (i.e.. 

m  P 

^m 

P      =   P.     ,    where      i    =   m    ,    defined   in   C      ) 
mi  p 

*m 


24 


V  -  Velocity   magnitude    of    the    fluid   particle 

V  -  Fluid  velocity   vector 
a    ,    <J>   -  Yaw,    Pitch   angles 

a.     ,    <J> .    -  Yaw,    Pitch    angles    for   probe    setting   i 

a_.     ,    <j>_.  -   Yaw,    Pitch    anales    for    the   assumed    flow   direction 
Ri            Ri 

direction    relative    to    the  probe    setting 

p             -  fluid   density 


25 


Bibliography 

1.  Dunker,    R.    J.,    Strinning,    P.    E. ,    and  Weyer,    H.    B. , 
"Experimental   Study   of   the   Flow    Field  Within    a   Transonic 
Axial    Compressor   Rotor  by   Laser  Velocimetry   and   Comparison 
With    Trhough-Flow   Calculations",    ASME    Journal   of   Engineer- 
ing   for   Power,    Vol.    100,    pp.    279-286,    April    1978. 

2.  Shreeve,    R.    P.,    Simmons,    J.    M. ,    Winters,    K.    A.,    and  West, 
J.    C,    Jr.,    "Determination    of   Transonic    Compressor   Flow 
Field  by   Synchronized   Sampling   of   Stationary   Fast    Response 
Transducers",    Symposium   on  Non-Steady    Fluid    Dynamics,    ASME 
1978  Winter  Annual    Meeting,    San   Francisco,    Dec.     1978.       (To 
be   published   in  ASME    Journal    of   Fluids   Engineering.) 

3.  Shreeve,    R.    P.,    McGuire,    A.    G. ,    and  Hammer,    J.    A.,    "Calibra- 
tion  of    a   Two   Probe   Synchronized  Sampling   Technique    for 
Measuring   Flows   Behind   Rotors" ,    paper   to   be   presented   at 
IEEE,    Eighth   International    Congress    in    Instrumentation    in 
Aerospace    Simulation   Facilities,    Naval   Postgraduate   School, 
Monterey,    September   24-26,    1979.      Published   as    IEEE    ICIASF 
Record   of  Proceedings. 

4.  Thompkins,    W.    T. ,    Jr.,    and   Kerrebrock,    J.    L. ,    "Exit   Flow 
From  a   Transonic   Compressor   Rotor" ,    AGARD   Conference    Pro- 
ceedings  No.    177,    Unsteady   Phenomena   in   Turbomachinery , 
pp.    6-1    to   6-2  3.      Meeting  held   at   the  Naval    Postgraduate 
School,    Monterey,    California,    22-26    September  1975. 

5.  Adler,    D.    and  Shreeve    R. ,    "A   General    Procedure    for  Obtaining 
Velocity   Vector    from  A   System  of  High    Response    Impact 
Pressure   Probes",    Naval    Postgraduate  School    Technical    Report 
NPS67-69-007,    July    1979. 


26 


°pendix  I  -  Program  VELOCITY 


VFLCC17Y  AND  DIRECTION 

;r    a    r:Luic    pj    fi    pgint   ugi.jg   fcur    pressures 


AUL      TAYLOR 


31 
PR 

AN 


iivrN    the    ;-<-'-: 

T  I  N  C  T     n  I  F  r.  C  T  I  i 

Ct;Ef'3    am:    flu  io    fi.,:pei5ti:s 

C     TUTAI.     AN:      STATIC     PRESSURE     AKF     CALCULATf 


su^cs  sfnjed  ey  FRceES  in  fcur 

N£t      AND     KNOWING     THE     ChABACT  ilU 'iT  ICS     CF      THi" 

THC     FLUID     VELOCITY.     DIRECTION 

T  A !_      AN:      fTAir 


)  I  v  'N  ■;  I  C  N 
"JCAL     nach 


PRfJl  (  19,2)  ,  FPFJ2(  19,2  )  ,  CP1  (  19  .  1 

AL  '  (  4  )  ,  PHI  (  4 )  .  NPRO  (4).PRESS(4) 


C  P  2  (  19 
(  4  ) 


I '"3  ) 


SUEKCUTINC      l^-UT      Fr.OEtVES     Th-:     NECESSARY     PPfl3r     Ct-ARACT 

rCF     THI"      T'.vC     i'PC-T.G      --     VATH  ICES     °  P  J  1     ANO     PPE2     RECEIVE 

THE     ax[S     (AL't-^     H,:     =>HI)     V>LUFS.     and     CP1      *NC     CP?     RECEIVE 

The  cp  ;aluc:  ron  probes  i  and  2  respect  ively  . 

CALL  INPT(N«LPF1  ,  NPHI  1  1  PRB1  ,CF1  , IER) 
..ALL  INPT<NALPH2.NPHI2.PRU2.CF2.IEP) 
IP  (  1  :  ■.  .NF  .  0  )  GCTG  10^0 


I  ST  I  CS 


WM 

GA 
TC 

cu 

50 

;0  10 


«LAC  f  F  E  F  C  L  L  C  V  I  N.  G  FLUID  FFCEFTICS: 

=  VCLECLLAR  HEIGHT  CI"  THE  FLUID 

N^A   -  FATIC  CF  SPECIFIC  HEATS  CF  THE  FLUID 

=  TEVFEK^TLFE  DEGREES  CENTICRACC  CT  Tl-E  CAS 

N-J    =  ESTIMATE  fF   TFE  CU  "PRE  S  S  I  C  I L  I  TY  FACTOR 


7  roo 


«  C  A  r  ( 

FCRVA 
PGA  £ 

*  r.  I  TE 

FOR  V  A 


1 


)00 


77  I  0 
C 

c 

c   ** 

c 


*P  I  TE 

ir  C  P.  N  A 


*  p  i  t  •-: 

F  C  fi  N  A 


j.'jcin 

KAF1C 

-  £3  14 
<  7  , 7  7  c 

T  (  •  FL 

•  PA 

'  :e 
( ^ .  e  c  0 
T(  •  1 

«  VE 

fX," 

(7,7  7  1 

T(///  . 

•    T  Y 


)  * 
.4  ) 
./W 
C  ) 
LID 
T  IC 

c  c 

C) 


N  ,  ",  AWVA.TC.CCVP 


hM,GAI»MA.TC.CCMF 
PROPERTIES  :»,//.•   yOLECULAP  wT  =  •  .  7  30  ,F  8  .4  ,  /  , 
CF  SPECIFIC  HEATS  -  •  ,  T  30  ,  F  f.  .  4  ,  /  ,  •  TEMPER  A  TUR  E  •  , 
'•iT33.rc.4f/.'  CCVPPCSS IE  I LTY  FAC7CF  3«,T3C.F8.4) 


i.nc 

C  ) 

•  p  f  c  n 
p 


STATIC         TC7AL,.12X.'VA'»      PITCH',  5  >  , 
ITY     MACH»  ,/,3X, •PRESS  (PA)     PRESS(PA>«, 
LE     ANGLE'.eX.Mf'/SCC)    KCMSEH*.//) 


J  Z  T  T  I  N  ! 


YAW        PITCH        PRESSURE'./, 
SETTING       READ  (PA)*,//) 


X*  *   START  LrCP  *  **  *>* 


RE  *D     IN      THE     EXPERIMENTAL     DATA    FOR     THIS     DE T f P V  I  N A T IC N 


NP 
AL 

PH 
PR 

AL 
AN 
TY 
TR 

TR 

10 

12  0 


15 

7720 
2  0 


Fe  (I)  =  Tt-F.     PRC5E     TYPE     (EITHER     1     Cfi     2)     CF     FFCCE     SETTING 

F(I)  =  ALPt-A      (YAW)     ANGLE     OF     PROBE     SETTING      I 

1(1)  =  FFI      (PITCH)      ANGLE     CF     FFCRF     SETTING      I 

£-5(1)  -  FFESSUFF     READ     fl  Y     PRC3E     SET7ING     I 


FI.      t:HJI,     NO^ni         CCNTATN     THE     N  S  Vi     VALUES     CF     THE     YAW,     PI7CH, 
C     PRCGE     7YPE     FCF     EACH     SETTING.      IF     THE     VALUE     READ     FCR     THE     PRC3C 
P^      I  _-     2LFL      (NPREII=0),     THEN     THE     PRCHE     SETTING     FCF     7HC     PREVIOUS 
•AL      13     USCC.        NC     DEFAULT     VALUES     ARE     PRCVIDEC,     SC     THE     FIRST 
IAL      »'<JST     CLNTAIN      THE     PROBE     SETTINGS. 


DC 
-SEA 
r  UP 
I  F  < 
ALP 
PH  I 
nll'il 
*RI 
FCP 

c:n 

Afi  I 


7730 


iO     I 

C(     Z     . 

i»A  T  ( 

(S  p  ',";  H 

(I   )  = 
(  I   )  = 

r(  I) 

T^  (  7 

V  A  T  ( 

T  I  NU 

TE  (7 

,v  4  T  ( 


=  1.4 

5C2  0 
3F  10 
I  .SO 
AL  F  I 
PI-  I  I 

=  nfi; 

.7  72 
1  >  .  I 
E 

,773 
//  ) 


,ENC  =  '59  9)     PRESS  (  I  )  .ALP  I  ,PH  I  I  ,  N^R'3  I 

.4,11) 

.0  )     CETC     15 


;'i 

1)     NPi-jrH  I)  ,ALP(I),PHI(  I).PRESS(I) 
j  ,?>  ,3^1 0.2  ,F1 4.2 ) 


C  ) 


ESTMLI3I-     SC'NNINC     RANGE,     GRIO 


WIDTH, 

27 


AND  INI7IALIZE  ERPCR 


COOOOO 1 0 

CC0000  2  3 

C  COOT  0.30 

C003040 

0  0  0  0  0^0 

zococe.z 
:  o  o  i  o  r  ■) 

cccocao 

C  0  C  3  0  "J  0 

r JOCO  100 

C  0  0  C  0  1  1  0 

"  C  0  C  0  1  2  C 

:  c  o  c  o  i  3  o 
c :  t co  140 
:coc j  150 
c  c  o :  o  i  t  o 

(*  C  3  0  0  1  7  7 

C  C  0  C  C  1  a  0 

:  :o  co  no 
:  oo  C0200 
:  o  o  c  o  ?  i  o 

:  c  o  o  o  2  ?  o 

C  C  0  C  0  2  3  0 
:  0  0  0  0  2  4  0 

C30C0250 

:  c  o  a  o  2  6  o 

0  0OC027C 
:  31072^0 
:C  3002Q0 
J0O0O  ?30 
"  C  0  C  C  3  1  C 
C30C0  320 
3  0OCO33O 
:30C034C 
C00C0  350 
C0OC0  360 
C3000370 
300CO380 
03  0  00  390 
3  0CC0403 
333  004  10 
300C04  2  0 
CCOC0430 
COO  00440 
CC0C3450 
C00C0460 
C0000470 
C00C0490 
CC0CC4T0 

coocosoo 

C0300510 
CCOC0520 
30300530 
C0000540 
C0OCOS5C 
00000560 
C00C0570 
CC000571 
37  0  00572 
C0OCO573 
COOC0574 
C3 3 00575 
C0300576 
03000580 
C30C0590 
C0000595 
00000596 
CC0C0597 
C0CC059e 
30000539 

coocoeoo 

0300061 0 
CC300620 
0  00C064C 
00  0  00650 
C3000660 
C0000C70 
30000680 


VELOCITY 


C         »1|:\,      JWJX     =      I'lMKC.      MAXIVuy     YA»     ANGLtG  J3  3C0C9C 

C        PM  lis  ,      T-VAX     -     vimkm,      VAXINUM     PITCH     ANGLES  00000730 

C        L3E  L                            =GFIC*ICTh  2  33 007  10 

C        ERFMIN                   =      ^IMM.iVl     ERROR     FCUNC     SC     FAR  33303720 

C  10000730 

AVU:-/,C.  C03C0740 

4VA>  =  .^;.  0  0030750 

P"U  =  -40.  3  C  3  C  0  7  fS  0 

PVA>=4  0.  0  0  0  C  0  7  7  6 

CE4_  =  5.  '01007H3 

130     tHT N  I  N  =  1C C CCO  .  10333790 

C  OOCCOdOO 

C                            ETAm     EC'N     FFCCECUfiE  C0OCC81C 

C  0  3300820 

CX=YA*     ANCLE     CUESS  C  0  0  C  C  8  .J  0 

C        Y      :     PIKh     ANCLC     CIFIG  C3320940 

C  s3  3OCO850 

Y=PMN  ^ICO^O 

1  t  0     X  -  A  v  I  is  C0  3Q0S73 

■ '  0  3  C  0  5  e  0 

C        CPS'JM        -     STl^Ec     THE     SUM     CT     THF     FClf;     CP«3     BEAD     FF.G*     '2  Y  CPC4L                                    '  C  3  C  0  8 -J  3 

C        PRSjUW     -     5TCFES     THE     SUM     CF     THC     FCUK     INPUT     PFESSUFES  33JC0900 

C        CPMK        -     STCFES     TeC     MINIMUK     CP     VALUE     FJR     THIS     GUESS  C0OC051O 

C         PKMN         =     STCFES     THE     PRESSURE     CORRESPONDING     TC     T  H  C     MINIMUM     CP  ^0000915 

C  C30C0920 

170CPSLV=C.  0  0  0  C  0  9  3  0 

P  R  3  S  U  V  -  c .  C  3  3  C  0  9  '»  0 

CPMIN-5.  0OC0SS0 

C  0  00  0  0  9  6  0 

C         ST/KT      THE     ANALYSIS     HY     TINOlNG     THE     (.f     VALUES     FSCV     T  t-  Z     T4ELC     CPCAL.  1 0  3  0  0  9  7  0 

C         A,MC     C/z-'LLATING     CPSU*.     CPMIN,     AND     PRSSUM  00030930 

00  3C09JO 

L3CZ0CK  =  1,4  3  0  0  010  0  0 

XR=>-ALF(K)  '3  0  0  C  1  0  1  0 

Y^=  Y-FH  I  ( K  )  CCO01 020 

IF(NPPE(K),E3.1)     CALL     CPCAKNALPhl  .NPhlt  ,PRC  I  ,CP  1  .XP  ,YR  .CP(K)  ,  IFL)  CCICI  3  30 

IF(  NPFC  <  K  )  .ra.2  )     CALL     CPCAL(NAL?H2  i  MPFI2  ,OR,;!2  ,CF2  ,XR  ,  Yfi  ,CP(  K)  ,  IFLJ0  1  )3  1  3  V3 

IK    lFL.NF.C)     CCTC     2H0  0  0  0  0  l  0  5  0 

CF3LM=CFSLM4CF<K )  0  C  0  C  l  C  6  0 

PRSSUM=PRS£LM*FRE<;S(K)  300C1 070 

IF(CJ-VIN.LT.CF(K))C'JTC200  3  C  0  C  l  C  fl  0 

CPVIN-CF(K)  3  0001090 

FF.N!N-PFESS(K)  3  0  0  C  I  1  0  0 

200     V.CNTINUE  00301110 

C  030C1 120 

FKC.V  TFE  AECVE  DATA,  CALCULATE  A  TOTAL  ANO  STATIC  FFESSUFC  00001130 

c  :  o  o  o  i  i  a  o 

C        {>[  i     -     fi     CHARACTERISTIC     TOTAL     PRESSURE     FOR     T  (- I  S     YAV,. FITCH  33301150 

C        PS  S     -      A     cHAF  ACTITMST  IC     STATIC     PRESSURE     H  OR     THIS     YAK,  PITCH  00301160 

c  :coc i 170 

3HI\CM=CFSLM-4.*CFiVIIN  0  30C  1  laO 

PTT  =  <     PRS SUM* <  1  .-CP^IN)      -     PRM  I  N* (  % .-CP S L W J      )/CENCM  303C1190 

,JL3-(      CPELV*PRNIN     -     PRS5UM*CPN  IN) /     DENCM  3C0C1200 

C  OOOCl Tl 0 

c   calcjlat::  a  characteristic  error  and  cjmfare  *ith  the  00001220 

C   PREVIGLSLY  rCLNO  SMALLEST  ERROR  C0301225 

c  00001230 

IF(  FTT.LE  .F5S)  CCTC  250  3  CO  01  24  0 

ERf  F=C.  00001  253 

J  G  i 2  £  I R  = 1  • 4  003012  6  0 

225  ERRF  =  ERFr  ♦  AfcC(CP(IR)  -  < PF E S S <  I R > -  PS S  > / < FT T-PS S >  )                  00001270 

lrr  f=.ej-rf  /n  .  coooi2ho 

IHtFMi.Cc.EHKUM     GCTO     250  C0001290 

C  30301300 

C    THIS  FCINT  HAS   ri-f  SMALLEST  CRHCR  HCUNO  SC  FAR,  SC   IT  IS  SAVEC          00301310 

ANT  FCFLACES  THE  PREVIOUSLY  FOUND  BEST  POINT  00301320 

C  CCOC  1  330 

C   PS  ,  PT  =  THE  EEST  STATIC.  TOTAL  FRESSURE  FCuNC  00001340 

XM  Its  ,   Y.MN  -  TFl   YAW,  PITCH  ANCLES  WHERE  THE  MINIMUM  ERROR  fcAS  FCUNC  C03C1350 

C  C0001 360 

LlihMN  =  EFFfi  C00C1370 

PS=FSS  00001 2a.O 

c  T  =  F  T  T  CC0C1290 

X  M  I  N  =  >  333014  3  0 

YMI  N  =  Y  00001 410 

'_  00001420 

250     X=XHDEL  C 0001430 

If  O.LE.AVAX)     CfTC     170  0  30  01440 

28 


VELOCITY 


300 

Y  =Y 
IF( 

c 

L 

At 

v-CN 

V- 

AN  i 

CLc 

L 

i  F  ( 

c 

l_ 

i*t 

9cP 

c 

LX  i 

LiiFT 

c 

AK  I 

^A 
I'M  I 
PM 
JEL 
GCT 

c 

c 

CAI 

ALL 

c 

1FL 

.j  TA 


4CEL 

>.LE.FM>)      ClTL      16G 

T1NCF      FEECCING     TPE     GPIO     SliE     oNTIl     TFE     E  P  P- t  P      IN     T  •"  E 
(-L  ACF.E  E      C.L     CEGRECS 

eel  .le  .c.'cd    cere    :ro 

••M      TFE      FfifCCCof-L     ARCLMJ     TK        Ej  E  S  T     F-CINT     F  C  L  IN  C     b'C     FAF 
US  INC     A     CHIC     1/3     AG     ..  «DE 

N  =  XK  IN-CEL 
>->f'Ih      *     DLL 
N  =  YVI  N     -     OF  I 
>  =  YW  IN      +     CEL 

=   n  c  l  /  2  . 

C      ISC 

tit     TFt     CFEIPFD     CJAMITIESi     FIRST     Ct-fCKINC     TCP     these     ep^civs 
*     o     MAhS     TPE     KANCC     CF      ThL     C A L I 3R AT  I  CN      TALLE     kAS     EXCSFDL'J 
1 f      TPE     LAST      SCAN 
"!IC     PfFESUPE     <  =     0,     TFE     FLCIC     VELOCITY     PECLIfiES     A     PCS1TIVF 
ST  AT  IC     FPESSUPE 


P.F  C      =  FLCTC  EENSITY         (K(J/M**2> 

VE  I     =  FLLIC  \ELCCITY     it'/ SCO 

CC         =  SCNIC  \TLCC1TY     CT     FLLIC      (V/SEC) 

M  A  C  h  =  FLLIC  MCI-     N  C  V  fc!  L  P 


,313  0  IFMFL.NE.C)  V-FITC(6,7CC0) 
7  0  CC  rLMAT('-M«  h/RNINC  THE  RANGE  CF 
I  MfELt  HChT  NET  FAVE  CEEN 

'/LLC»  FFCPEf;  CAICCL AT  ICNS  •  ) 
IF  (  C£  ,LE  .C  .  )  CCTO  450 
PEC  -  P£/(FCAEHLVP*nC^72.1f.)) 

n,acf  =  sgpt(<  (pt/p:j*  +  (  (gamma-i.j/cai»ma)-i 

c  0     =     SCF.T(G^NMPGAi*<TCt272.1t)J 

vLL  =C  C*VACF 

*P  I  IE  (C .6  C  IC  )     FStPT.XVIN.YyiN.VEL.NACF 

FChfAK  IX  ,*F12.2.5X,2Fe.2tSX,Fe.2.F3.3) 

G  C  T  C      1  C 


THE     C  ALI  E  F.  A  T  It  N 
SJfi    IC   IENT      TC      • 


)/  (  (  GAMA-  1  .  )/2  .  )  ) 


CO  10 

C 

t    A 

L 

•'+3  0 
70  10 


CCCATIVE  STATIC  FPESSUPE  FAS  8EEN  FCUNC 


1 


1CCC 
70  JO 


99 


*P4 
FCP 

CCT 

k.  P  I 
FEK 

SIC 
ENC 


1 


TE(£.7C1C)  T£  ,PT  ,XVI  N  ,  YMIN 

MT(i   NEGATIVE  STATIC  PHF.SSJKE  •  ,/  i 

•  P£  ,PT  ,  YAfc  ,PI  TCF   :«.4F12.2> 
L   10 

I L  ( t  .  7  C  -  C  ) 
NAIC  /f.   INFLT  CRPCP  CCCCPPEO  *HILC  •  . 

"PEAEINC   IN  THE  PPCCE  C  H  AR  AC  T  E  R  I  S  T  IC  S  •  ) 
F 


cooo 
cccc 

CCOO 

cccc 

CC  JO 
CC  JC 

cooo 

c  co  c 
c  coo 

CCOG 

cooo 
ceo  c 

C  00  c 

ccoc 

co :  c 

CC  J  0 

ccoc 
ccoc 

00  JO 

cccc 
ccoc 
cooo 

CO  JO 

CCCO 
COOO 
COOO 
CCCC 
CCOO 
CCOC 
CCOC 
COOO 
CCOC 
COOO 
COOO 
COOC 
OOJC 
COOO 
CCOC 
COOC 
CCOC 
COOC 
COOO 
CCOC 
COOC 
C  0  JO 

cooo 

COOC 
00  00 
CCOO 
JCOO 
0000 
CCOO 

cooo 


1  450 
146C 
1470 
1  480 
1  490 
150  0 
I  SI  0 
1  520 
I  S30 
154  0 
1  550 
1  560 
I  570 
1  5^0 
1  590 

1  coo 
i  e  10 

1620 
16  3  0 
1  6*0 
16  45 
1650 
1660 
K7C 
1690 
I  700 
1710 
I  720 
1  730 
1  740 
1  750 
1  7o0 
1770 
1  780 
1  790 
1  30C 
1810 
1  820 
1830 
1  840 
1  850 
1  860 
1  900 
1  91  0 
1  920 
1920 
1  940 
1  950 
1  960 
1970 
1  980 
1  990 
2C  10 
2020 


29 


VELOCITY 


C  :C0C2C3C 

C  C0O02040 

JJJl- .:i.  1  [  KC   [NPT  00002050 

c  :coc2C6o 

C        Tt-  :li      UJr'XlINr     l!E/>DS      It.'     TFE     CATA     fTi;     THE     PRCEE     Ct-  A  3  ACTE  R  I  3  T  I C  3 .  JC002Q70 

C        IT     k.A\     HI"     Ch/NGLD     TC     ANOTHER     3L  I  T  AiJl_  IZ     FORM      IF  RIICL'IREO.                                               J00020d0 

c  J  CO  C2  C30 

<-      »*i  -:    =  nuitep.  cr  points  on  the  al-ha,  phi  axis  :ooo2i)o 

^                    "II.  (>,i)       -      ALPHA     VALUE3    CN     ThP     A  X  I  "5     nr     THE  <-  «  L  I  F  R  A  T  I C  N                                           nDC'lM 

<~                                                        T  A  C  L  E  G  0  0  C  2  1  2  0 

C                  lJKJ(*i2)      =     PHI      VALUES     CN     THE      AXIS     CT     THE     CALIOPATICN  00002130 

c                                     taele:  J0002140 

C                  CP(  NA  ,NP  )  =     KATPIX     CCNIAIMNG     TFF     VALUES     CF  CF     FCP     THE                                             00002150 

C                                                       PARTICULAR     PPODF  50002160 

C  0  0  0  0  2  17  0 

C  C  0  0  0  2  1  3  0 

SLJFCLTINF.   IFPT  (  NA  ,  NP  .TkJ  ,CF  ,  IER)  2  5  0  C  ?  1  9  0 

:iVENCICh  FFHII9.2),  CP(19,19)  G  3  G  C  2  Z  5  C 

KCAC(t'.eCC0.FNC  =  S9S>NA.NP  C03C2213 

dJJO  KL.;NAT("H)  CC0C2220 

i>i "^(>  .ecic.'.Nc-oi'))    (p^Ki.n.i-i.NA)  :000223a 

REAC(i:igC10.rNC  =  1?5)      (PP3(  J,2)  ,  J=l  ,NP)  0  0  0  0  2  2  4  0 

(20  10     l-CRMKlCFStS)  J  0  0  C  2  2'j  G 

RE.  A  C  (  t:  .  9C  SO  ,LNC'  =  <}SS  )      (  (  CP  <  I   ,  J  )  .  J  -  1  .  .>IP  )  .  I  -  1  .  N  A  )  -10  1  0  2  260 

2C20     FCR  VAT  (  1  CFE  .5  )  3  0  0  C  2  2  7  0 

I E  R  =  C  28302280 

:<LTLKN  :01C?29  0 

i.  0  0  3  0  2  J  :  0 

Tt    II.  H     FA£      fl  FN      AN     C^PCK     WHILE      INPUTTINC     T  F  E     CATA,  3  0  0  C  2  2  1  C 

C                     3C     AN     Fl- r  O     FLAGi      IER.      IS     SET     =1  G0002320 

C  :O0C22  1C 

•GQ5  IFU=1  0  CO  02340 

RETURN  0  00  0  2  350 

L.\D  5  00  02  36  0 


30 


VELOCITY 


<-  : :  v  c  2  3  7  o 

C  :C3C2390 

C            SUUP-CLTINH  CPCAL  )30  32  3<S0 

C  :33C2430 

C  NA  ,NF  -  1  CF  MPH<  ANC  Oh[  ANGLES  IN  THE  CF  C*LtCP*TICf,  '.'  3  0  C  2  4  I  C 
C  Alt*)  -  VALUES  GT  THE  ALPHAS  FCR  THE  C  AL  I  t3R  A  T  I  C  N  TAOLE  (YA.v  ANCLES)  OOJ02420 
C   P<rP)   -  VALUES  OF  Tl-J  PHIS  ECP   THE  CALIOPATICN  T«L;LE   (PIICI-  ANCLES)    "3032430 

C   CP(NA,r>P)  =  VALUE  CF  rp  FOR  EACH  ANGLE  SET   (  MNA).P(NP)   )  J300?'*40 

C    X        "  ^LSIFfC  AlEHA  ANCLII  <".  0  3  0  2  A  5  0 

C   Y        -  CESIF-EC  PHI  ANGLE  23302460 

C   <!       =  CALCILMFC  CP  0  0002470 

I   IF  LAG  =  KRFCK  FLAC  3  00  02  48  0 

v.  C  3  0  C  2  4  ;  0 

C   T'li:  rFCGSAN  EST  I  NATES  THE  VALUE  CF  CP  Kn  4  C  I  V  L"  N  ANGULAR  I\.'-T  2 3 0  02 50  0 

C  (ALPHA,  rt-  I  )  USING  A  LINEAR  OC'JPLE  I  NT  E  H  3  GL  A  T  I  C  N  SChcVF.  F?fT*EEN  ;  0  0  0  ~"  Z  1  0 
v.    TF  E  KNC'aN  VALLFS  CF  CP  FCT  ANGLES  At!CVE  ANC  EELCU  THE  C  '  S  I  n  E  C  A\CLC    23002:523 

C  )  0  0  0  2  5  J  0 

c  ;o-)or50 

S  U  a  F  C  U  T  I  N  i"  C  P  C  /  L  (  N  A  ,  N  F  ,  F  R  13  .  C  ^  .  X  .  Y  .  ?  ,  I  FL  A  C  )  '.'  0  3  0  2  5  5  0 

i3IVENSICNFFil(NA,2),CP(NA.NF)  13332  5  60 

l.  0  C  0  0  '  5  »  3 

C   GT*«r  THE  SEAFCH  FCR  THE  ALPHA  VALUES  A U C V E  ANC  CCLC*  T  F  C  20302583 

G   OEElRLiC  YA'a  ANCLE  0J102590 

:  o  o  o  o  2 1 )  o 

vxa,  NNA  -  STORES  THE  ENTRIES  TO  P  I- S  AN"  CP  FCF  THE  ANCLES     033C2C1C 

G                     lUCVr,  Jfc'LCW  ThE  3PSIRF0  ANGLE  0  03  02  620 

C              \r>  ,   ak  =  THE  ALPHA  ANGLES  ADCVEi   CCLCa  Tt-E  3ESIRCC  ANGLE  C0302fi:!0 

C  2C0C2640 

J  C   1  0  I  -  2   i  N  A  0  3  0  0  2  6  5  0 

v  X  A  -  I  0  0  3  0  2  6  6  0 

M  N  A  =  I  -  1  2  0  0  0  2  €  7  0 

A--Fr;C(I.I)  0000258  0 

*.\=FRii  (  N'NA  ,  1  )  0  33C269  0 

I  'T  (  ' n  .  S  G  .  X  .  A  NC  .  A  N  .  L  E  .  X  )   G C  T  C  2  2  30002700 

10  :CNIINUP  00  3  02710 

,_  30002720 

G         IE     Thu     LfCP     t-A^     CEFN     COMPLETED     WITHOUT     FINCINC     ANGLES     SLPRI3UND  ING  C0002730 

G         THE     IOCS  I  RFC     ANCLE,     THEN     AN     ERRCR     FLAG     --      IFLAG     --      IS     SET:      IFLAG=1  30002740 

C  30002750 

IFLAG=1  COOC2760 

RE  F  LKN  3  30  C2  7  70 

C  0  0  3  0  2  7  e  c 

C    XO  -  F  PAC  T  ICr- Al.  CISTANCE  CF  THE  CESIRED  ANCLE  2ETUEEN  THE  300027^0 

C         KNCaN  C*LI!.?FMICN  ANGLES.  00002800 

0  0  0  0  2  S  1  0 

25  XE=  (  X -AN  )  / <  AP- AN  )  00002820 

c  00002830 

C   THE  SC/RCH  FCR  THE  PHI  VALUES  STARTS.   VARIABLES  A  F  E   IDENTICAL  C0OC2840 

G    TC  THCCE   IN  THE  PFFVICUS  SEARCH  EXCEPT  »P«   SUBSTITUTES  TCR  'A«.  C0002853 

C    ANC  »Y»  ANC   'J*    FFJLACE  •X*  AND  M*  PCSPECTIVFLY  00002855 

.-  C30C2860 

JC   ::  J=2.NF  00032870 

VXP-J  00002180 

VNF=J-1  C0OC2893 

rr-rrflj'l  0G0C2900 

p^hrlmvnfU)  22221212 

IF<  FP  .GE  .  Y  .  ANC.PN.LE.Y)  GOTO  45  00002520 

.-r-»  ■.  T  ,.  ■  .r-  COOC2<330 


'- 


30     CCMIM'E  2222oc^2 

1FliS-  =  1  00002*540 

RETLRN  0C0C2550 

45     YJ=(Y-r>N)/<FP-FN)  2222o^2 

00002Q70 

C         *i:     NG'*     FINE     THc     V/LUES      IN     THE     CP     CALIL5RATICN     T3ELE     AHICH     CCf.RESPCNC         CC0C2S80 
TO      THt      (_ALiLHATirN      ANGLES     AECVE     AM     3ELCW     THE     JESIFED     YAW     AND     PITCH         00002<5<iO 

00003000 

Cll=CP(VNA.l*NP)  orin'lo^o 

C12=GP(VN/.v>r)  2 

C2  ">  =Ch(  V>  A  ,  »»XP  )  w0003C-+0 

CC003050 
C0003060 
3C003070 
00003C80 

Z  =  XE*YC*<  C22  +  C  1  1-C12-C2  1  )  +  XE*(C21-CU)  ♦  YH*(C12-CM)  +C11         300030O0 

CC3C3  100 
C    IN  IERP CLATICN  SUCCESSFULLY  CEVPLFTED,  ERKCF  FLAG    IFLAC=0  00003110 

V  U  ' J  V  J  I  t  V 

IP,  «r -r  C0003 130 

^V.ch  300  03  140 

g^  .  .  COOC3150 


c^  1  -zr  (  vx  t , VN0  ) 

Thc   1NTHPPCLA7LD  CP  VALUE  eETWEEN  THE  TCUK  KNCHN  CP 
VALUCS!  Cll.  CI  2.  C2  1.  C22 


G         V ALUCS 

c 


c 


APPENDIX  II 
VELOCITY  NOTATION  SUMMARY  -  main  program 

ALP (I)  -  Yaw  angle  of  probe  setting  I 

AMIN ,  AMAX  -  define  the  minimum  and  maximum  yaw  (alpha)  angles 
of  the  search  grid 

COMP  -  the  compressibility  factor  of  the  fluid 

CP(K)  -  C   interpolated  from  the  appropriate  calibration  table 
for  p  probe  setting  K 

CPMIN  -  stores  the  minimum  C   found  during  this  guess 

P 

CPSUM  -    stores    the   sum  of   the    four   C    's    read  by   Subroutine   CPCAL 

P 

CPl,  CP2,(I,J)  -  C  calibration  table  for  probes  1  and 

C$  -    sonic  velocity 

DEL  -  search  grid  spacing  (degrees  of  angle) 

DENOM  -  stores  an  intermediary  mathematical  quantity 

ERRMIN  -  stores  the  minimum  error  found  so  far  for  the  problem 

ERRR  -  C   average  error  characteristic  for  the  guess 

GAMMA  -  ratio  of  specific  heats  for  the  fluid 

IER  -  input  error  flag  =  0  means  no  error,  =  1  an  error  occurred 

while  readinq  in  the  C   calibrations 
^         P 

IFL  -  interpolation  error  flag  =  0  interpolation  accomplished 

=  1  range  of  the  calibration  table 
was  insufficient 

MACH  -  fluid  M-ach  number 

NALPH1 ,  NALPH2  -  number  of  yaw  angles  across  the  edge  of  the  C_, , 
Cp~  calibration  tables 

NPHI1,  NPHI2  -  number  of  pitch  angles  across  the  edge  of  the  C  -.  , 

Cp2  calibration  tables 

NPRB(I)  -  probe  type  for  probe  setting  I  (either  1  or  2) 

PHI (I)  -  pitch  angle  of  probe  setting  I 
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PMIN,  PMAX  -  define  the  minimum  and  maximum  pitch  (phi)  angles 
of  the  search  grid. 

PRBl,  PRB2  (N,J)  -  contains  the  alpha  and  phi  angles  for  use  with 
Cp,,  C  ~  respectively.   J  =  1  refers  to  yaw  angles 

J  =  2  refers  to  pitch  angles 

PRESS (I)  -  pressure  read  at  setting  I 

PRMIN  -  stores  the  pressure  at  the  setting  corresponding  to  CPMIN 

PRSSUM  -  stores  the  sum  of  the  four  input  pressures 

PSS  -  contains  a  static  pressure  characteristic  for  this  guess 

PTT  -  contains  a  total  pressure  characteristic  for  this  guess 

RGAS  -  ideal  gas  constant  (Joules/kg-  K) 

RHO  -  fluid  density 

TC  -  fluid  temperature   C 

VEL  -  fluid  velocity 

WM  -  molecular  weight  of  the  fluid 

X,Y  -  yaw,  pitch  angle  guess  (one  of  the  search  grid  points) 

XMIN,  YMIN  -  yaw,  pitch  angle  where  the  smallest  error  was  found 

XR,  YR  -  yaw,  pitch  angles  of  the  guess  relatwe  to  the  probe 
setting  being  considered. 
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NOTATION    SUMMARY    -    SUBROUTINE    CPCAL 


AP,AN    -    Yaw   angles    above    and  below    the   desired   yaw   angle 

CP(NA,NP)    -    C      calibration    table 
P 

Cll,    C12,    C21,    C22    -    C      values    surrounding   the    desired   C 

P  P 

IFLAG  -  error  flag  =  0  means  the  interpolation  succeeded 

1  the  range  of  the  C   table  was  too  small 
y  P 

MNA,  MNP  -  Stores  the  location  of  the  calibration  yaw  (alpha) , 
pitch  (phi)  angles  below  the  desired  yaw  and  pitch  angles. 

MXA,  MXP  -  Stores  the  location  of  the  calibration  yaw,  pitch 
angles  above  the  desired  yaw  and  pitch  angles. 

NA,  NP  -  number  of  yaw,  pitch  angles  in  the  C   calibration  table 

PP,  PN  -  Pitch  angles  above  and  below  the  desired  pitch  angle 

PRB(N,K)  -  Contains  the  yaw  and  pitch  angles  for  the  calibration 

table 

X,Y  -  Yaw  and  pitch  angles  where  a  C   is  sought 

XB,  YB  -  Fractional  distance  of  the  desired  yaw,  pitch  angle 
between  the  known  calibration  angles 

Z  -  the  interpolated  C   value  for  X,  Y 
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NOTATION  SUMMARY  -  SUBROUTINE  INPT 


CP(I,J)  -  Calibration  table  read  from  the  file 

NA  -  Number  of  yaw  angles  on  the  edge  of  the  C   table 

NP  -  Number  of  pitch  angles  on  the  edge  of  the  C   table 

PRB(N,K)  -  contains  the  yaw  and  pitch  angles  for  the  C   calibration 
table  p 

K=l  yaw  angles 
K=2  pitch  angles 


35 


APPENDIX  HI  - 


-91.00 

-'10  .00 

-  7  C.  CO 

-'•0  .CO 

-  '2  i .  .  0  0 

-40.00 

-  30  .  C   3 

-  ?  o  .  :  3 

-icon 

0.0 

1  o .  o  o 

.'  c .  ;  o 

J  C  .  0  0 

4  0  .  C  0 

c  ).n 

(    3.30 

70.00 

«  3  .  3  0 

i  C  .  2  ? 

-  O.J) 

-  j  J  .  0  J 

-  7  C.  CC 

- 1  0  •  c  c 

-5  COO 

—  -4   0  .    J  >J 

-  i  ■:, .  o  > 

-2  0.      1 

-13.   >  0 

1  .0 

10.  00 

3  3.30 

3  C  .  JO 

4  0  .  C  0 

5  0.00 

<:,  0.30 

70.  OO 

►«  0  , 

S  0  #   . 

).  CIO  J? 

o .  o ;  o  o  j 

0  .   )  C  0  C  0 

3  .  )  0  C  C  3 

0  .  C"   1  3  3 

)  .    3  0  0  0  0 

0  .  0  )  J  0  3 

0  .       30  3  0 

3.303     C 

0  .  .3  J  2.   3  0 

1 .000  JO 

0. 00  00 J 

o .  o  co  oo 

0  .  J  0  c  c  c 

0  .  C  i  j  0  0 

J  •  0  C  C  C  0 

3  .  3  0  J  3  0 

o  .•■■•:-  3 

3.333     0 

3  .  .03 0  )  0 

J  .  C  0  0  «  3 

0 .  10  359 

0  .  0  C  7  <_   ; 

0  .  0  1  ?  (  r 

0  .  C  1  7   15 

0.0  27  7  9 

0 .026  '2 

0,0  2  .  3  ? 

3  .  0  7j3'J 

'  .   3  2-433 

3  .  JOo  7  2 

0  .  0  2  2  7  j 

0  .  0  1 7  d  'j 

0  .  C  1  2  f  C 

O.CO7'     . 

0 .  ;  C  35  ) 

0 . 3  2  0'     i 

0  .  3  0  3  0  C 

v    •      .     »-                  v 

0  .  )0  2<    I 

0.01  4  26 

o  .  o  2 :  j  _ 

0.04*77 

0  .  0  7  C  1  9 

.3  .  C  V  )\  1 

3  .  1  C  -4  .» t: 

3.114"1 

'.3.1  1  7     2 

0  .    1    14  2  1 

o.i  )  -*  7 :, 

3  .  0  3-u  I 

0  .  0  7  C  l    ) 

0  .0  4s  7  7 

j  .  o  2  c : .' 

3  .  C  1   4   •>>■ 

)  .  2  0  2  6  < . 

3 . 3 3  0  30 

o.)c  : , c 

1 .  :  "  j  i  3 

o .  r>  J  1  7  e 

J.  0'»(     t  3 

) .  t  :  j  j  j 

0.15227 

0.1  >.  2f  '. 

0.22'*J5 

0.24557 

3,25271 

C  ,24  ":-  7 

0  .2  24  ■   '. 

0  .  1  9  2  <  ' 

0.1  >2  J  7 

3.1CS.J 

0.06   40 

0 .02 1 2- 

3  .  C  0-11  3 

0  .  0  C  C  0  0 

3  .  0  0  0  0  0 

0.01   3=3 

3  ,  C  5  3  2  1 

2.1:2  :<o 

).   1M.'7C 

0  .  2  I  5  0  J 

0.3217, 

J  .  3  7  4  4  C 

0  .  4  l.  9  3  H 

0.   4  1'"    5 

0 ,40 3  0* 

C  .  77  14  0 

0.3217  ' 

3.2   i552 

J  .  1  S  J  7 9 

J  .  1  1  2  2  ) 

0.  0  52  21 

J  .  3  12  6   ! 

O.CCOJu 

0  .  0  C  C  ,3  u 

0  .  0  2  u  0  -J 

o.c  '  n  l 

3  .  1  t  2  2  J 

3  .  . "    ''  "-    i 

)  .    3 '    O  t    5 

3.4CC32 

0  .  '  j  24  2  0 

0  .'-;ti  1  2  4 

..  .  .  '    .'  ,0 

0  .  5  rt  1    '4 

0  .  2   '  4  2  ' 

0  .  4  t  C  7  2 

3  .  3  t  ^  f  >• 

■3  .  v>  1  2  . 

3 .  i  <  2  ;•;  -i 

0.07711 

J  •  0  ^  i.  0  -i 

>  .  C.  '  C  c  0 

3  .  C  3000 

0.    '2  >95 

3.0095    ! 

0.23  i.37 

3  .  3  2  •'  4  ..' 

3,47   -  7   ' 

o  .uj  .. .; » 

o  * '.  ^  2  ;  j 

j  .  7  4  _'  2  'j 

)  .  ,' '    2  7  7 

0  .7  4,.'   i9 

3.i    220  .3 

0  .  5    3  9  3  4 

0.4703  9 

3  .  !  3  8  4  2 

3.72  .07 

J  .  o-.  9  5  S 

o .  o  2  j  ;  o 

0  .  0  C  J  u  0 

0  .  10CCC 

J  .  C  1  0  4  9 

3.1   K. '    0 

0  .  >  a  '  .  4  ' 

">  .  3  <    7  14 

1  .  J  5  1  ^  i 

.3  .  6  -i  O  •)  j 

0.810.3  3 

0  .    IC^Jj 

J  •  d  5  J  c  t„ 

0..4  69  6  ': 

0  .  fc  C  3  . !  3 

0.6  9  0  J.S 

0  .5   -17  1 

3  .   3 f    7  1   '4 

3.2454* 

3  .  1  1  (  w  0 

J  .0  50  4'/ 

o .  o  :  o  j  o 

0  •  )  C  0  0  0 

)  .  0  J  H.1 

0.  1  "■':  M 

0  .2c(   3  3 

3.43142 

0  .  <-  2  3   )( 

3  .  752  3H 

3,37215 

3  .  5  4  84  j 

0  .    !  7  4  i  7 

0  .  -^   4  J  4  J 

C  ...  72  lf" 

3.7'  ."   ••• 

0  .6  0 COo 

3  .  '.    3  1  <4  2 

/      „      1      i 

2  .  1  2  '  •  "»  I 

3.33-30? 

).G)  J>) 

) .  j  :  c  ;  j 

v..  o  ;~  _ , 

0  .  1  2  e  2  1 

3.2^027 

3  .h  3  v;  ! 

0.*    1 .    11 

) .  r  t  ,  j  j 

3  .  ,934  ? 

)  .9 7  2  78 

1.00000 

0  .   )  ?  2  7  J 

0  .  6  ?  ;  .4  7 

o .  tc;"' 

J.cl.Ml 

'.'  .  v2*i9  7 

0  .    ."J1' 

2  .  1  2  •=  J  1 

0  .  0  '     '   7  9 

3 . 0 2  3  3  0 

>  .  C  3  0  0  0 

0.  0  2  1  /J 

0  .  1  22  20 

0  .  ,">  c3  c  4  2 

0  .  4  2  330 

0.59)   T4 

0,746   3 • . 

0  .8(    '"  1 

3.94753 

0,9  '4  57 

J.    M7V1 

0.86  5:31 

0.7-4     CC 

0  .   i  9  I  7  4 

0.422   !0 

0 . 3  5  942 

0.1    'l^J 

0.231  ;  3 

3.33   300 

3.30)30 

J  .  0  *85c 

0.11022 

0.22  :  u  4 

0  .  3H222 

>  .  r  ;•  7  j,o 

0  .   >»01  7 

3.7<-.<4  <  .4 

0  .  ■-'  (   ^   ;  4 

C . d 9  3c' 

0.  3-.  624 

).  7  3454 

J.OEC17 

0.5  2  73d 

0  .  3  n  2  2  2 

0.22  364 

3.1 10'3 

0  .  02.156 

0.3033 3 

J  .  3  3  3  3  0 

3 .02   >7U 

3  .  )  y,!00 

o  .  i  -;  j  7  j 

0.3  2 1 39 

0  .4  54 C  1 

0 .57^9" 

0 .67624 

3  .  7<*06  i 

0  .  7  e    J  7  / 

3.74053 

0  .  5  7i  24 

3  .  5  7  u  is) 

0  .  4  L  4  J  I 

0.3  21  j  9 

0.1957 C 

0  .  C  -   -1  C  0 

3.02!  .3  7  9 

0  .  3  3  C  3  0 

). 3  1822 

3. Q 7  06  1 

0.  1    636  . 

o  .  2  -t  -3  j  o 

0.3?.     29 

3  .  i  «4  <  4  C 

0  .52-31 

0  .-v  791.8 

0,597  5)  0 

0.5  7  0--  3 

3.52831 

0  ,44'J4C) 

0  .    !  3.:  2  v 

0  •  2  4  £3  0  > 

o.i  5 j t : 

3.07 C t I 

0  .  C  1  *  2  2 

o.ocooo 

J  .  '3  0  3  3  0 

0.0122 7 

0.  C  4:3  80 

0 . 1 0439 

J  .  1  /  J  7  1 

1.     ''43.3 

O.J  1 4  <.'  1 

0  .  3  7  C  2  7 

J.lC'.  )2 

3 . i 196  5 

)  .40C  0  2 

0.37027 

3.31421 

0.24  ">  5  8 

J.  1  7.17  I 

0  .    1    J  4  0    , 

0.  C41    3  0 

0  .  <:'  l  ?  5  7 

C  .  C  C  -.  3  0 

0  .  3  C  C  0  0 

3  .  0  0  7  4') 

0.0291 3 

0.3'   334 

0.10313 

0 .  14  7  2  1 

J  .  1  J  3  j  9 

0.  2  22  c '3 

0  .,'44  Cf 

0  .  2  '..?-'  1 

0  .  2  4  4  ?  < 

0 . ?T 268 

0  .  1  t  Ht.9 

3.14721 

.3  .  1  0  "■  3  3 

J  .    362  3  4 

0.  02S  1   3 

0  .   )  C  '  4  , 

3  .  3  0  C  C  0 

3  •  C  C  0  0  0 

0  .  0  3  2  4  0 

3.01351 

3.0  2  D) 

0.04801 

3  ,0684  1 

3.0.)/76 

0.  1  2  3'.  1 

0.1    14  0  0 

0  .  1  1  7  f  z 

0.1140 0 

0.10261 

0  .  0  °  7  7  6 

0 . 0  f  3  4  3 

0.04301 

3  ,0209  5 

0.01 351 

).i)03rj 

0  .  3  C  0  3  3 

0  .  toe  CC 

0 . C 0  0 no 

0  .0 3  34 H 

0  .  0  0  7  4  6 

0.0  1 2  3  7 

3,017f  I 

3.072  5 0 

3.  0  2(  6 2 

0  .  0  2  <5  2  -3 

3,0202: 

3  •  0  2  c.  2 .« 

0  .  C  2  f  r>  2 

0.02  2  5  ', 

3.01761 

3.21237 

)  .  A  C  7  '4  t 

0  .00  34  8 

J  .  0  0  0  9  0 

0  .  0  J  0  J  J 

0.0  CO  J  J 

O.OOCC 0 

3.CO0OO 

J  .  3  C  0  3  0 

3  .  C  3  0  3  0 

J. 3  3330 

C  ,  0  0  '3  0  0 

3.300OO 

0.  COO  0  J 

o .  o :  j  3  o 

0  .  0  C  0  J  0 

0  .IOC  CC 

0  .  C  0  3  3  0 

3.  3C00O 

0.0)3  30 

0.30330 

3  .  3  C  3  3  0 

3.30000 

3.CJO0J 
10         10 

-  0  )  .  0  0 

-  Hi  0  ,  JO 

-  7  COO 

-r;0  .  CO 

-5'    .0(3 

-4  0.00 

-.3  0  .no 

-3  0  .  3  0 

-13.00 

0.  0 

1  J  .  J  0 

2C  .  J  J 

2  C.  20 

4C  .CO 

5  0.0  0 

o  0  •  0  0 

7  0.  3  0 

8C.3Q 

0  0.00 

-90 .00 

-30  .00 

—    *  C  •  <J  sy 

-t  C.C3 

-50  .   30 

-43.00 

-30.CO 

—  ~»  P          *  c"1 

-  1  C  .  0  0 

0.0 

10.00 

2C  .  3  0 

3  C  .  C  0 

'4C.CC 

5  0.00 

^C.OC 

7  0.00 

ci  0  .  3  3 

'50.00 

0 . C  0 0  JO 

0. 0  3  00  J 

0  .  0  C  0  0  0 

0  .  0  C  C  C  0 

0  .  C  c  0  0  c 

o.  :  c  Cuo 

O.C  3  J  3  0 

0  .  O  J  3  ■"  Q 

0  .  0  0  3  0  0 

0.03)30 

0.0)000 

0 .  0  0 J  0  0 

o .o  c  ico 

J.OOCCJ 

0  .  C  C  0  3  0 

0  .  J  c  0  0  0 

0  .  C  C  0  0  0 

0.3  3  33  0 

0  .  0  C  3  0  0 

3  .  30  0  30 

J.O^O'l  J 

0  .    12  2  5  0 

0.  J  C  7 1.  2 

0  .~12£G 

2 . C  1  7  35 

0.02  2  7  5 

3*02672 

0.22'-:  *  C 

3  .  0  2  3  2  C 

2 .3  2^30 

0  .  0  '  6  7  2 

0  .  0  2  2  7  5 

u  .  0  1  7  a  '2 

3  .  0  1  ?  f  C 

0  .  C  C  1 6  5 

0  .  3  C  3  5  9 

c  .  n  3  C  0  3 

0  .  3  0  3  3  0 

3  .  C  C  0  0  0 

0  .'3  0  36S* 

J. 01 *26 

0 .  J  J  0J2 

0  .  0  4  0  7  7 

0  .070  19 

0  .  0  "  9  1  1 

3 . 1 C43S 

0.11421 

0.11762 

3.11*21 

3 . 1 0  4  35 

0.0-391  1 

0. 07C  1  ) 

0.0  497 7 

C  .  C  2  C  2  2 

0 . C  14  2  6 

0 . 0C3l 9 

3.03030 

3.30300 

3  .  3  C  4  1  3 

0.03135 

0 .066  4  0 

0  .  1  C  .3  3  J 

0.1  =  .10  7 

0  .  1  9  2  c  5 

0  .  22  4  35 

0.24557 

3.25271 

0.24557 

0.22485 

0  .  1  9  36  6 

0.  1  32  3  7 

0  .  1  C  £  -3  2 

u.O  £<   4  0 

0.-33  1  35 

0  .  C  0  3  1  3 

0  .  0  C  0  0  0 

0  .  C  0  0  0  0 

0.01382 

3.05221 

0.11210 

J. 19279 

0.  2=5  3.' 

0.32179 

0  .  3  7  4  4  C 

0. 4C20H 

0.  *  196  5 

0. 406  08 

0.37*40 

0.72  1 79 

0  .2  3552 

j.  rjc'7'j 

0.  1  1  2.;  J 

C  .0  5  !2  I 

0.01362 

0  .  C  0  2  0  0 

0  .  3  C  0  C  0 

0.0200H 

0.0  7711 

C . 16220 

0  .36  3  25 

0  .  3l  t.6fl 

0  .  4  c  C  2  2 

0.52423 

0.58124 

0  .  5  9  75  0 

0  .  5  0  1  3  4 

0  .534  2  3 

0  .  4 1  3  2  2 

0.26668 

0.2t  326 

0  .  16228 

0.07711 

CO  20  C8 

3  .  0  C  C  C  C 

C  .  C  C  C  0  0 

0. 02595 

3 • 0  •  •  95  j 

0.2  0907 

3  .  3  3  ii  4  3 

0 .470  39 

0.  2-3  3  34 

0 . t    3  2/0 

0.74239 

0  ,7t.  27  7 

0  .74  2  39 

o  .  c  e  2  o  o 

2  .5  8  9  7<4 

0  .  4  7  C  :  9 

3  .  73t,4  3 

0  .  30-40  7 

0.  0')  ?5  J 

0  .  0  2  5  9 :5 

0  •  C  C  J  C  G 

C  .  0  :  C  C  C 

C. C 7049 

0  .  1   16  9  0 

0.24  54  6 

0.39714 

0.-5171 

0  .69093 

3.0303  3 

0.3t5!j 

o .a  <5J6t 

0  .  fi  6  C  £  £ 

C  .  i  C  0  J  3 

0.69093 

0.55  1 71 

0.3':  7  1  <4 

3.24549 

0. 1 1 o90 

0.0  3049 

O.OOC)  0 

0  .  J  C  0  J  0 

0  .  0  3  2  C  2 

0  .  12<,  M 

0  •  ?  e  6 .33 

0.43142 

0.6033c 

0.7522c 

3.87215 

).  '5484  3 

0.  9  7  4 '5  7 

0.9  43  4  3 

C.H72  15 

0.75  2  2  0 

0  •  6  C  C  0  6 

0.43142 

0  .  2  '  ■  f  3  2 

0 . 1 26  7  1 

0.03303 

0 . OOOOO 

0  .   3  0  3  J  J 

U.O  2o 20 

0  .  1  2  8  2  1 

0.27027 

0.42893 

O.C  1?  1  1 

0  .  7  f  9  2  3 

0.39247 

0.97278 

l.oomo 

0. 97  2  ? 3 

0 .8  5347 

0  .  7  6  ?  2  3 

3  .£  121  1 

0 .4  2  8 93 

0.27  0  2  7 

0  .  12831 

0  .02239 

0, 30030 

0  .  C  0  0  0  0 

0.02  100 

0  .  1  £  2  d  0 

0  .25? 4  2 

0  .  4  2  2  (3  0 

0.39174 

0.  7^6  0r. 

0.86331 

Q.9*75C 

0.97467 

0.9  47  50 

o .  m  e  a  a  i 

0  .  7  1'jOo 

0  .'^  1  7  4 

0.422^0 

0 .2594  3 

0.12280 

0  .  0  2  1  < ".:  C 

C.CCCOO 

0.00 1 0  0 

0.02 356 

J  .  1  1  C  2  2 

0.2  236  4 

0  .  2  P  2  2  2 

0  .  5  7  738 

0  .f  HOI  7 

0  .79454 

0.8f  -324 

C  .  8  9  3  €  6 

0.96824 

0 . 79454 

0  .  6  e  0  1  7 

0.32722 

C  .  2  '3  2  2  2 

0.22  364 

0.1  1022' 

0.02856 

C  .  30030 

0  .  C  C  3  0  0 

0.02379 

0.0  i  ?  0  0 

o .  i <5  5  7  o 

0.3  2  139 

0  .  4  5  4  C  1 

0  .5  759 Q 

0  .67634 

0.74053 

0.7C  2T  7 

3.74053 

0.67624 

0. 57699 

0.45  401 

0.321'39 

C.  1  957C 

C .Cc200 

0 . 32379 

0.03000 

0.30330 

0 .0  1  322 

C  .07061 

0.1    3  J  ':  t 

0.24333 

0.3  5.  2  9 

0.44c.tC 

0.52331 

0.3796 0 

3.59  750 

0.  57*568 

3.52831 

0  .  4494  0 

3.  2  5229 

0.2463 5 

0. 1 5ucc 

0  .  0  7  C  t  1 

0  .  C  1  6  2  2 

0  .  0  C  0  0  0 

3.30000 

0.01357 

0 . C4330 

0.  134  39 

0  .  17  2  7  1 

0.24  :U 

-.31421 

0  .  37C2"7 

0  .1G(    3  2 

0.4  1 965 

3.406^2 

0  .37027 

0. 3  142  1 

0  .24553 

0  .  17  2  7  1 

0  .  1  3  4  3  9 

0.  0  4.    6C 

0.0  1257 

0  .CCC  30 

0  .  3  C  0  0  0 

0  .  0  C  7  4  9 

0.0251  0 

0.06  234 

3.13333 

0 .  14  7  2  1 

o .  i  s  3  c  o 

J  .  2  2  _  C  r? 

0  .24  4 Cf 

0.15271 

0 .2449f 

0 .22268 

0.13369 

0. 14  721 

0 . 10  313 

0  .  01234 

0.02^  10 

0  .  0  C  7  4  0 

0  .00 CC 0 

0  .  C  C  C  0  0 

0  .  3  0  3  4  fi 

0. 3  135  1 

0 .02595 

0 . C4801 

0. 06843 

3.0<    776 

0.  13  3,1 

J  .  1   14  0  0 

0  .  1  1  7  e  2 

0.  1  1400 

0. 1 C36  1 

0  .  0  £  7  7  6 

0  .  C  fc  3  4  3 

3  .  C  4  3  0  I 

0.02395 

0  .  0  1  3  '.5  1 

0 . )JJ4J 

0  .0  Co  00 

ij  • «)  C  w  c  c 

0  .  C  C  C  9  0 

0  .  0  C  3  4  6 

0.30746 

0  .  U  1  2  3  7 

2. 0  I  761 

0.02256 

0.0  26  62 

0.  02 O^J 

0  .  0  2  0  2  0 

0  .02<-:  2  2 

0  .  C2f.  b2 

0  .0225C 

0  .  2  1  7  e  1 

0.01227 

3  .CC  746 

3.00348 

0  .  C  0  0  9  0 

0. occoo 

o.  J  c  ;  CO 

0  ♦  JCCCO 

0  .  C  C  C  3  0 

0. OCCOO 

3.  cooco 

0  .  C  0  3  3  C 

O.CC  300 

0. OOOOO 

0.010  0  0 

3.0  CO  )0 

0  .  0  C  i  i  0  0 

0  •  ,)  0  C  C  0 

0  .  C  C  0  )  0 

0  .  0  C  C  0  0 

C . 000  30 

0  .  0  C  3  3  0 

0.30300 

0.03000 

0  .  C  0  0  0  0 

^■~ 

SfWPLE  INPUT  (f.;on*-n 


28  .80 
1 0  7660 
.i  0589  0 

j  o ;;- !-;  ^  ii 

•!.  I.)  7 ': .  8  0 

1  ■')  9  !.  4  a 
10  9740 
J.  0  7  n  9  it 
1  i.  39  »Q 
10PBP.0 
I  095  9() 
i  0  9  !i  90 
105980 

981  0  0 
1 0545(1 
1  0  ;:i  '  2  0 

995  6  0 
1.200  00 
i.  1. 5  0  9  0 
1  1  c>  0  1  0 

i. :  lv;  n  j  n 
1 l 8840 
j 1 0640 
1 16650 
11 1370 


i  .40 

-PS  ,0 

20,  p 

n  . 

l.Q 

0  . 
2%  .  0 

ii  , 
0  . 

., 

u . 

as ,  o 

•".> 
*... 

o . 

0  . 

0 . 

P  5 

.  o 

-;-'S .  o 

0  . 

2  5  .  0 

0  . 

37 


APPENDIX  *JZ  -Sample  OuTpuf 


fluiij    fkopcrtil: 


M  C  L  il  C  U  L  A  K     AT      - 
;?a  r  I  )    cf    -,pl  C  [ I-  [c 
TtVPhlR  ATUttC     rcc     C     = 
CCMPRr.GLl  lILTy     FACT  L 


(■  l  \  i  :    = 


.•f'.rfCCC 
I  .  4  C  C  0 

:  c  .  k  <:  o 
I  .ec  r o 


•RL.ir 

T  YHK 


i  T  T   I  .\  j 


!J  I  TCH 
S  H T  I  I  N  G 


r:(  at    ( p a  ) 


-  2  £  .  C  J 

<-.  .3 

2  5  .00 

c  .  c 


L  .    C 

c . : 
c .  c 

J  .  CO 


1  C  7  C  C :  j  .  0  0 

icce<:c.oo 

I  C  7  5  £  C  I  0  C 


■2  5  .0  J 

C  .0 
2  £  .  C  2 

C  .0 


0. 

0  . 

r  . 


CO 


1  C  0  1  <  0  ,  C  0 

ico7*c.co 

1  C  7  0  ?  C.  C  0 

i  eas c :  .co 


l 

-  2  Z   .  'J  J 

i. .  ( 

I  C  2  1  I  C  .  C  0 

I 

C  .3 

c .  r 

!  C  '?  1  c  c  .  C  0 

1 

25  . :  j 

c.  c 

ic<;  i  q  c .  cc 

M 

C  .0 

.2  : .  C  0 

1C59€C  .00 

I 

-25  .00 

<:  d  l  c  c .  c  o 

1 

C  .0 

IC5450.G0 

1 

.2  5  .  C  0 

c .  c 

IC4  120.CC 

p 

0.0 

,:  5  •  c  j 

c  c>  i  £  c  . :  o 

I 

C  .0 

c  .  C 

I  2  0  C  C  C  .  0  0 

2 

C  .0 

2:.C0 

1  1  5  0  'CQC 

1 

-25  .J  J 

C  .  C 

115C1C.CG 

1 

2  5  .00 

0  .  c 

1  15  0  1  C  .  C  0 

L.J 

C  .0 


u  .    C 

2  5.C0 

C  .  c 
0  .  C 


i  i  a  e  a  e  .  c  o 

1  1  0  (5  4  C  •  0  0 

i  i  e  e  e  : .  o  c 

1  1  1  2  1  0  •  3  C 


pni 


.STATIC 
0  3     (  P  A  ) 


icr^ 

^;;C3  5  (PA  ) 


/NCLE 


r^  ITCH 

ancl: 


VLLC^ITV     VACH 
(V/SEC)    NCVCCh 


i  000  t<)  .00 
30  0  19.10 

9  I  0  C  4  .  1  0 
eS£  79.94 


1  I  G  1 5  i.  •  ?  5 
HCOSf.l? 
ll0l2V.bc 
1  10  1  J  C. 25 
1  i ,9£ £.  7 j 
I  2  0  1 3  €  .  2  c 


:c.;7 

— ?  .  1  5 

1  2  .  f  0 

2  2.15 

-C.CC 
-  7  .  C  4 


J  0  .  J  7 

5  .  C  £ 

0.^7 

-  le.52 

—  O.JO 

-  ^  .  <..  2 


I  2  d  •  7  6 

1  2  fi  •  j  J 
187.4J 

1  £  Q  •  0  3 
220. o  7 

2  2  6.  ? o 


.  J74| 

•  57? 
.54  4 

•  ~  4  5 
.  C  4  1 
.f  c  7 


38 


APPENDIX  V 
NOTES  ON  THE  USE  OF  VELOCITY 

INPUT:   The  required  input  consists  of  probe  calibration  data, 
fluid  properties,  and  finally  the  experimental  pressures.   Sub- 
routine INPUT  reads  the  calibration  data  from  each  probe  type 
in  the  following  form: 

1.  The  first  card  contains  the  number  of  yaw  and  pitch  angles 
on  the  axes  of  the  calibration  table  (format  214) 

Ex:   19   19  means  19  yaw  and  19  pitch  angles  were  used  in 

the  calibration  and  the   C    table  will  therefore  be 

P 

19  x  19  in  size. 

2.  The  next  few  cards  contain  the  values  of  the  yaw  angles 
where  calibration  points  were  taken  in  the   C   table. 
Values  are  entered  in  format  F8.2,  one  angle  every  8  columns 
After  all  the  yaw  angles  have  been  read,  the  pitch  angles 
are  entered  starting  on  a  new  card. 

3.  The  experimentally  determined   C  's   of  the  calibration 

surface  can  now  be  read  for  each  angle  pair  starting  from 

the  smallest  yaw  and  pitch  angle  and  with  the  pitch  angle 

varying  most  rapidly.   Ex. 

C  (-90).  -90).  C  (-90,  -80)...  C  's  are  read  format  F  8.5. 
p      '         P  P 

All  of  the  calibration  data  are  read  on  Machine  Unit  8: 
Cards  are  assumed  to  be  80  characters  in  length. 

The  following  fluid  properties  are  entered  next: 

Molecular  Weight 

Ratio  of  Specific  Heats 
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Fluid  Temperature  Deg  C 

Compressibility  Factor 
Machine  Unit  6  reads  this  data  from  one  card,  Format  4  F10A. 

At  last  the  experimental  results  are  entered.   Four  cards 
are  required  for  each  trial,  one  card  per  setting.   For  format 


Columns   1-10 

11-20 

21-30 

31 


Experimental  pressure 

Yaw  Angle 

Pitch  Angle 

Probe  Type  (1,2,  or  blank) 


If  Column  31  is  left  blank,  only  the  experimentally  read 
pressure  is  registered;  yaw  and  pitch  angles  for  that  setting 
remain  unchanged  from  the  previous  trial.  The  first  trial  must 
contain  angle  settings  and  probe  type  since  no  default  values 
have  been  assumed.  Again  machine  unit  6  is  used  to  read  this 
data.  When  no  more  experimental  pressure  data  is  available, 
the  program  terminates. 

The  experimental  pressures  can  be  based  in  any  absolute 
system  of  measurement;  ex.:   Psia,  KPa,  Atm,  mmHg,  with  the 
same  numerical  results  (the  units  in  the  titles  of  the  static 
and  total  pressure  columns  will  not  apply) .   The  analysis  below 
shows  that  in  determining  velocity,  the  pressure  units  cancel. 
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The  velocity  is  calculated  from   V  =  MCq  ,  where,  from  Eq .  (1), 


M 


Y-l 


(PT/PS) 


-1 


1/2 


(Y-D/2 


and 


Here , 


Cq  =  /yRT 


C   =  sonic  velocity 


M  =  Mach  number 
Ps'PT  =  fluid  static,  total  pressure 
R  =  ideal  gas  constant 

8314  Joules/kg  mole°K 
MW 

T  =  Fluid  Temperature   K 

V  =  Fluid  Velocity  (m/sec) 
p  =  Fluid  density 

Y  =  ratio  of  specific  heats 

CPCAL:   A  linear,  double-interpolation  scheme  is  employed  to 
determine  a  value  of  C   between  four  points.   A  second-order, 

IT 

double-interpolation  scheme  has  also  been  devised  and  tested, 
and  is  presented  at  the  end  of  this  report.   Figure  V-l  is  a 
graph  of  the  accuracy  of  both  schemes  as  a  function  of  the  num- 
ber of  calibration  points  in  the   C    table.   Values  were 
determined  by  filling  a  calibration  table,  extending  from 
-90°  to  +  90°  in  yaw  and  pitch  with  the  C  's  which  would  result 
from  an  ideal  probe,  and  testing  6084  points  (78  x  78)  within 
the  table.   If  no  highly  unusual  distortions  in  the  calibrations 
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of  the  probes  occurs,  Figure  V-l  shows  that  a  significant  reduc- 
tion in  the  amount  of  calibration  required  is  possible  with  a 

second  order  scheme.   Further,  if  the  accuracy  of  the   C 

P 

determinations  is  known,  Figure  V-l  can  provide  an  estimate 
of  the  number  of  points  needed. 
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Second  -Order     Double- Interpolation    Scheme. 


C  0CCC447C 

C  0CCC44eC 

c  r^cc44<;c 

SLERCUTINE    CPCAL  (  SA.hF.FRB.CP.  X,  Y,  2.  IFLAG  )  OCff^OC 

C                            NA.NP     =     *     OF     ALPHA     AND     PHI      ANGLES     IK     THE     CP     C  AL  I  PR  A"r  IFN  0  C  C  C  4  5  1  C 

C  A(NA)      =     VALUFS     OF     THE     ALPHAS      IN     THF     CALIFSQATICN     T  »BI    F      (YAW     AN<"|   orcr4r?C 

C  P(NP)      =     VALUFS     CF     THE     PHIS      IN     THE     CALIBRATION     TABLF      (PITCH     AN  CI   0  C  C  C  4  5  3  C 

C                            CP(NA.NP)     =     VALUE     CF     CP     FCR     EACH     ANGLE     SET      (      A(NA),^>(NP)      )  0CCC4540 

C                             X                  =     DESIRFD     ALPHA     ANCLE  CCCC45CC 

C                           Y                 =    DESTRFO     FH     ANGLE  CCCC45fC 

C                             Z                   =     CALCULATED     CF  0CCC457C 

C                            IFLAG     =    ERROR     FLAG  0CCC4560 

C  THIS     PROGRAM     ESTIMATES     THE     VALUE     CF     CP     FOR     A     GIVFN     ANGULAR      I  NP<   0  C  C  C  45  <;  C 

C  ALPHA.     PHI.     USING     A    LINEAR    OCUELE     INTERPOLATION     SCHFWF     BETWEEN     K NCOO C C 4 C C 0 

C                 VALUES    OF     CP    FOR     ANGLES     ABCVE     ANC     EEIC*     THF    DESIRED     ANGIF.  OCCC4eiO 

DIMENSION     PRE(NA.2>      .CP(NA.NP)  0CCC4f?0 

DO     10     I=2.NA  1CCC4f?C 

20     MXA=I  0CCC4f40 

MNA=I-1  0CCC4CeC 

MMA=MXA+l  OCC04ff«- 

IF  (MMA.GT.NA  )     MMA=VN/-1  ZCCC4t?e 

APaPRB(I.l)  0CCC4ffC 

AN=PRB(MNA,1 )  CCCC4f7C 

AO=PRB(MMA. 1  )  ncCC4f7c 

IF  (AP.GE.X.AND.AN.L E.X  )     GCTC     25  0CCC46EC 

10     CONTINUE  OrcC*ec:C 

IFLAG=1  0CCC47CC 

RETURN  O0CC471O 

2=5  DO  30  J  =  ?.NP  0CCT47?C 

40  MXP=J  0CCC474C 

MNP=J-t  0CCC4750 

MMP=MXPfl  0CCC475C 

IF (MMP.GT.NP )  MMF=MNF-1  0CCC47C« 

PP=PRB(J,2)  0CCC47CC 

PN=PRB(MNP,2 >  0CCC477C 

PQ="PB(MMP,2)  CCC0477C 

IF (PP  .GE. Y.AND.PN.LE .Y )  GCTC  45  0CCC47fC 

30  CONTIN'""  0C0C47«^C 

iflag=i  occr^ecc 

PFTUPN  0CCC4P1C 

45  Cl i=CP(MNA,MNP)  0CCC463C 

CI 2=CP(MNA.MXP)  CCCC4P4C 

Cl 2=CP(MNA,MMP)  0CCC4P45 

C21^CP(MXA.MNP)  OCCCAE^O 

C22=cp( mxa.mxp)  crrc4pfc 

C23-CPIMXA.MMP)  0CCC4e70 

C3  l-CP(MMA.VINP)  CCCC48E0 

C32=CP(MMA.MXP)  0CCC4e<-C 

C33=CP(  VMA.MMP)  0CCC4<;CC 

F1-(X-AP)*{X-A0)/(AN-AF)/(AN-AC)  OCCCAOIO 

F2=(  X-AN)*(X-AO)/( AP-AN  )/  (AP-AC  )  0CCC4<;?C 

F3  =  ( X-AN)*(X-AP)/( AC-AN)/( AO-AD)  0CCC4^3C 

CI =F  1*C1  1+F2*C21*F3*C31  0CCC494C 

C2-F  l*C12+Fr*C22+F3*C32  OCCC4<5cC 

C3=F1*C13*F2*C23>F3*C23  0',CC4<;«C0 

Z^(Y-PP)*(Y-PQ)/(PN-FF)/ (PN-PQ  )*C1     ♦  CCCCAS7C 

X         (  y-PNI)*  (Y-PO)/(P°-FN  )/(PP-P0  )*C2     ♦  CCCC4<5PC 

2         (Y-PN)*(Y-PP)/(PO-FN  )/(PC-PP)*C3  0CCC4<;c0 

RETURN  0CCC5C00 

END  CCCC501C 
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